GCM study of hexadecapole correlations in super deformed 9 Hg. 
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Abstract 

The role of hexadecapole correlations in the lowest superdeformed band of 
194 Hg is studied by self consistent mean field methods. The generator coordi- 
nate method with particle number projection has been applied using Hartree- 
Fock wave functions defined along three different hexadecapole paths. In 
all cases, the ground state is not significantly affected by hexadecapole cor- 
relations and the energies of the corresponding first excited hexadecapole 
vibrational states lie high in energy. The effect of rotation is investigated 
with the Skyrme-Hartree-Fock-Bogolyubov method and a zero range density- 
dependent pairing interaction. 
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I. INTRODUCTION 



Recent experimental data on the spectra of superdeformed (SD) bands in 149 Gd [|TJ, 
153 Dy 0], 194 Hg [§], and 131,132 Ce [§] have put into evidence a Al = 4 staggering of the 
dynamical moment of inertia of these bands. It manifests itself in systematic shifts of the 
energy levels which are alternately pushed down and up with respect to a purely rotational 
sequence. The amplitude of this staggering is of the order of 50 eV. Since these oscillations 
distinguish states differing by four units of angular momentum, it seems natural to explain 
their origin by a coupling between the rotational motion and hexadecapole vibrations []I]||. 
This scenario requires the presence in the mean field of the nucleus of a small component 

7T 

with a four- fold rotational symmetry, i.e., a symmetry associated with a — rotation around 

one of the principal axes of the quadrupole tensor. The strongly prolate-deformed intrinsic 
Hamiltonian is then slightly perturbed by a term with the C4 symmetry. 

The effect of such terms has been investigated in phenomenological models with the C4 
symmetry axis coinciding either with the symmetry axis (z-axis) of the quadrupole tensor 
0-^] or with the C4 perturbation quantized along the rotation axis (x-axis) |]||. Hamamoto 
and Mottelson [|] have studied the properties of a quartic rotational Hamiltonian. Staggering 
appears then as a result of tunneling between the four equivalent minima of the total energy 



surface of the Hamiltonian due to the AT-mixing. Burzyhski et al. [|1(| have obtained the 
exact solutions for a system of iV identical particles in a single j-shell interacting through 
a multipole-multipole Hamiltonian with quadrupole and hexadecapole terms. They have 
shown that the staggering in 

J(2) 

may occur in certain cases. 
All these studies have assumed the presence of a hexadecapole term in the Hamiltonian 



of the SD nucleus. It still remains to explain the origin of such a term. Ragnarsson [11 
has considered triaxial hexadecapole deformations in a study at high spins of SD bands in 
the A=150 mass region using the shell correction approach with a Nilsson potential. The 
effect of £44 deformation has been found extremely small. Frauendorf et al. |12[ have applied 
the tilted axis cranking formalism also to a Nilsson Hamiltonian with a hexadecapole field, 
Y44 + 14-4, and have concluded that it is impossible to generate a staggering of the observed 
order of magnitude even assuming 644 values as large as 0.1. 

In the present work, based on self-consistent methods with Skyrme interactions, the im- 
portance of hexadecapole correlations in the lowest SD band of 194 Hg is investigated. Firstly, 
we analyze dynamical hexadecapole correlations without rotation by means of the generator 
coordinate method (GCM) with particle number projection |R|. Secondly, cranking calcula- 
tions have been performed within the Hartree-Fock-Bogoliubov + Lipkin Nogami (HFBLN) 
method [|I4|,|I5| with a zero-range density dependent interaction M] in the pairing channel. 



II. THE METHODS 
A. The Generator Coordinate method 



The GCM formalism used in this study has been presented in detail in Ref. [13]. It 
permits one to study the dynamics of the nucleus as a function of one or several collec- 
tive coordinates using generating functions with good particle number. The intrinsic wave 
functions, |$), are obtained by constrained HF+BCS calculations with the Routhian 
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E' = - A($|iV|$) - A 2 (AiV 2 ) - qUi^lQ^) - Q 44 ) 2 . (1) 

In Eq. (|l|) N is the particle number operator and Q\ 4 denotes the hexadecapole moment 
quantized along the i axis, i.e., 

Q\, = \j r\Yl, + YU)p(v)d\ (2) 

where i = x, z and p denotes the selfconsistenf HFBCS density distribution. The functions 
Y44 and F 44 are the spherical harmonics along the x and z axis respectively. Thus Q\ A 
has a fourfold (C4) symmetry along the zth axis. The energy dependence on hexadecapole 
distortions around the SD minimum of 194 Hg is determined using constraints on either Qf 4 
or Q 44 . We have chosen z as the symmetry axis of the quadrupole tensor, whereas x denotes 
the perpendicular axis. 



In the particle- hole channel, we employ the Skyrme interaction SkM* | 17|| . In the particle- 
particle (pairing) channel, the state-independent seniority interation is used. An approx- 
imate variation after projection on particle number is performed by means of the Lipkin- 
Nogami method 



In this study it is assumed that the density distribution has three symmetry planes. The 
HF equations have then been solved on a rectangular mesh 12 x 12 x 16 fm 3 , corresponding 
to the octant with positive values of the three cartesian coordinates. Single-particle wave 
functions are discretized with a mesh size of 1 fm. We have checked that the results are not 
qualitatively modified when increasing the size of the box or changing the mesh size. 

From the HFBCS wave functions, states \NZ) with the correct number of particles are 
obtained by projection: 

\NZ)=P NZ \$), (3) 

where 

Pnz = - f exp U(j) N (N - N)\ d<j> N - Texp {i<f> z (Z - Z)\d<\> z (4) 

TV J0 *- ' IX JO > 

is the projection operator on neutron and proton numbers. The integrations over the gauge 
angles are approximated by n-point trapeze formulae. The choice n = 7 ensures that the 
particle number dispersion never exceeds 0.001. 

The GCM wave function is a superposition of projected wave functions corresponding to 
different hexadecapole deformations: 

I$>gcm = / f(Qu)\NZ;Q l u )dQl 4 (5) 

The weight function / is obtained from the Hill- Wheeler equation: 

J [(NZ; Q\ A \H\NZ- Q 4 ' 4 > - E(NZ; Q\ 4 \NZ- Q 4 ' 4 )] /(Q 44 )dQ 4 ' 4 = (6) 

where H is the many-body Skyrme+pairing Hamiltonian. Finally, the collective wave func- 
tions, g, can be obtained by an integral transformation of /(Q 44 ): 

g(Q\d = J (NZ; QUNZ; f (7) 
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B. The cranked HFBLN method 



The HFBLN method for rotating nuclei has been presented in Refs. Mffa . As in the 
GCM calculations without rotation, the SkM* parametrization of the Skyrme force has been 
used in the particle-hole channel. In the pairing channel, the density-dependent pairing 
interaction, 



V P = ^(1 - P a )5(n - f 2 ) 



ap 



r\ + r 2 



(8) 



has been employed. The strength Vq has been chosen equal to 880 MeV and the reference 
density, p c = a -1 , to 0.16 fm -3 . It has been shown recently [IBJ that the introduction 
of a zero-range density dependent pairing interaction improves considerably the agreement 
between the HFBLN calculation and the experimental data for SD bands in the v4=190 
mass region. In particular, the experimental transition energies in the SD band of 194 Hg are 
reproduced with a systematic error lower than 10 keV over 12 transitions. 

Hexadecapole constraints have been introduced in the cranked HFBLN equations in the 
same way as in the nonrotating case. 



III. RESULTS 

The GCM collective path is determined by a full minimization of the energy with con- 
straints on the hexadecapole moment alone, as defined by Eq. (]!]). This does not guarantee 
that other collective variables do not change along the path, i.e., that the hexadecapole mode 
is decoupled from all the other modes. In the case of such a coupling, one should in principle 
enlarge the collective subspace and include other variables as generator coordinates. 

In our GCM calculations with the Q Z 4A constraint, the quadrupole- hexadecapole coupling 
is very weak. Namely, the mass quadrupole moment, Q2o=Q2o =r2 ^20' does not vary by 
more than lb in the whole considered range of Q Z AA values. This is not true for the Q^ A 
trajectory. Indeed, the Q% 4 moment can be expressed in terms of hexadecapole moments 
Ql^ quantized along the z-axis: 

Q44 = qQaa "I ~^~Qi2 "I -^t-Qao- (9) 

Consequently, Q44 has p=0 and 2 components along the z axis which are strongly coupled to 
the quadrupole moments Q?, and Q\ 2 - To avoid many- dimensional GCM calculations in the 
Q44 direction, we have defined two different trajectories. In the first version, the quadrupole 
moment varies self-consistently as a function of Q44. In the second version, the quadrupole 
moment is constrained to its value at the SD minimum. 

The potential energy curves corresponding to projected wave functions @ are plotted in 
Fig. [I] as functions of Q\ A (top) and Qf 4 (bottom). The evolution of the quadrupole moment 
Q20 as a function of Q44 is plotted on Fig. ||. In the same figures are shown the results of 
the HFBLN calculation performed at I =42. 

According to the calculations, the collective motions along the Q44 and Q\ A paths are 
decoupled. Namely, in the Q\ A version of the calculations, the value of Q\ A does not vary 
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by more than 1000 fm 4 along the collective path. Also in the Qf 4 version the value of Q44 
remains always compatible with zero. 

The potential energy curves shown in Fig. [I] exhibit similar dependences as functions of 
the hexadecapole constraints. This means that, qualitatively, the response of the nucleus to 
the hexadecapole perturbation of the average field is the same in both cases. As expected, 
relaxing the constraint on the quadrupole moment (solid line in Fig. |l|) leads to a slight 
decrease of the energy. In spite of the fact that the quadrupole moment shows large variation 
along the unconstrained path (see Fig. |2|), the corresponding gain in energy never exceeds a 
few hundred keV. 

The effect of rotation on total potential energies seems to be fairly weak, although the 
conditions of the calculations with or without rotation are not exactly the same. The curves 
are surprisingly close in the Ql 4 variant. In the calculations with a constraint along the 
rotation axis, there is a slight shift of the minimum toward smaller hexadecapole moments at 
7=42 but the hexadecapole stiffness is very similar. The small shift in the equilibrium value 
of Q44 with angular momentum can be attributed to the rotation-induced shape change in 
the SD yrast band of 194 Hg. As discussed in Ref. |19[, a systematic decrease in both a Q20 an d 
Ql with rotational frequency is expected for SD yrast band in 194 Hg. According to the self- 
consistent calculations, the expectation value of Q\ 2 is almost zero (i.e., at least two orders 
of magnitude smaller than Q z m ). Consequently, according to Eq @, the equilibrium value 
of should follow closely the value of Q\ . Indeed, the shift in Q z m from approximately 
38000 (1=0) to 35300 (1=42) fm 4 is consistent with the results displayed in Fig. |. 

It is worth noting that the quadrupole moments minimizing the energy for small values 
of Q44 are rather different for 7=0 and J=42 (see Fig. 0). This is probably more related 
to the different pairing interactions employed and to particle number projection than to 
rotation. Indeed, the quadrupole moments of the unprojected HFBCS states |$) with and 
without rotation vary in a similar manner as functions of Q44. On the other hand, the 
neutron pairing energy which is of the order of 4.0 MeV for both pairing interactions around 
the SD minimum, increases up to 6.5 MeV for a seniority interaction at small Q44 values. 
This leads to large differences between the projected and unprojected wave functions. 

A static calculation does not give access to a scale parameter which would permit one 
to compare directly the evolution of the energy as a function of Q44 and Q\ A . As discussed 
above, the large values of Q% A are related to the large hexadecapole moment of the nucleus 
along the z axis, while the mean value of Q44 at the SD minimum is compatible with 0. 
A first crude measure of the variation of the collective wave functions along the collective 
path is given by the variation of the overlap between the wave functions corresponding to 
different values of Q 44 . According to calculations, the overlaps between the wave functions 
corresponding to the SD minimum and the extreme hexadecapole deformations considered 
(7?*~4MeV) vary between 0.1 and 0.3. These rather large values imply the rigidity of SD 
states with respect to the hexadecapole fields with fx=4. 

A better way of comparing the collectivity along the collective paths considered is to 
perform the full GCM calculation. This calculation gives at the same time the distribution 
of the lowest wave function as a function of hexadecapole deformations and the excitation 
energies of the associated modes. We have performed GCM calculations along the three 
collective paths defined above. The Q% A path with a constrained quadrupole deformation 
has been used to test the discretization of the hexadecapole moments, with up to 17 wave 
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functions corresponding to hexadecapole moments ranging from 4000 to 36000 fm 4 . These 
tests have shown that eight to ten discretization points lead to an accuracy of the order of 
100 keV, which is sufficient for our purpose. 

The results for the three collective trajectories are summarized in Table |. The HFBCS 
SD state is predicted at an energy of 1519.41 MeV and its hexadecapole moments Ql A 
and Q44 are and 20000 fm 4 , respectively. These values are not strongly modified by 
the collective correlations. The predicted gain in energy due to dynamical correlations for 
the collective lowest state is in all cases lower than 1.0 MeV, which is much smaller than 
the typical energy gains associated with the quadrupole and octupole degrees of freedom 
(see for instance Ref. [0). Also the quadrupole and hexadecapole moments of the first 
excited collective states are not significantly different from the values obtained for the static 
minimum. The energies of the first excited state are large, i?*>3MeV. The slightly lower 
value obtained in the unconstrained Qf 4 path is probably more due to the variation of the 
quadrupole moment along the path than to a genuine hexadecapole effect. 

In Fig. || are shown the squared moduli of the collective wave functions g, Eq. (0), cor- 
responding to the lowest GCM state along the Q| 4 (top) and Q44 (bottom) trajectories 
discussed above. It is seen that the collective ground states do not show any significant 
distortion with respect to Q44. The distributions are centered around the static minimum. 
The collective wave function obtained in the unconstrained Q44 calculations is slightly asym- 
metric, i.e., components with lower hexadecapole moments have a slightly larger weight, but 
this effect is not large enough to modify the total hexadecapole moment. 

The GCM results that we have discussed show that hexadecapole correlations are rather 
weak. In particular, no hexadecapole excitations at low energies have been found. However, 
the wave functions considered do not belong to a irreducible representation. To construct 
such wave functions, one should consider the mixing of the four wave functions obtained by 
rotations of — around the quantization axis of the hexadecapole tensor. For instance, for 

the C4 symmetry associated with the z axis, the four corresponding transformations of the 
HFBCS wave function |$) can be expressed as permutations of coordinates: 



(x, y, z) => (y, -x, z) => (-x, -y, z) 



-y,x,z) 



(10) 



Let us recall that if one imposes the condition that the nuclear density is symmetric with 
respect to the three planes x = 0, y = 0, z = 0, one can decompose the single particle wave 
function into four components with definite parities with respect to these planes PT 
Namely, the wave function $^ can be written as: 



V $/M / 



im$ fc (f,+) 
Re$ k (r, -) 
V Im<f> k (r, -) J 



(11) 



This choice ensures that the parities of the components §k,a with respect to plane symmetries 
about x , y and z planes depend only on the label a and the parity p k as indicated in Table 
|l[ Up to an irrelevant phase factor, the permutation of x and y in wave function (|ll]) gives 
rise to (for positive signature): 
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( Re$ k (r,+) \ 

Re$ k (r,-) 
\Im$ k (r,-) J 



( - J Re$ fc (F,+) \ 
Im®k(r r ,+) 



-/m$z,^r ,, , 



V 



(12) 



0/ 



where r 1 is obtained from r by the permutation of x and y. Using relation (O) and Table II, 
one immediately sees that the second and fourth symmetry operations in Eq. (|10D lead to 
wave functions with spin-down component orthogonal to the spin-down component of 
provided that the $>k,a have axial symmetry with respect to the z-axis. The overlap between 
the original HFBCS state and the total wave function obtained by the third transformation 
of Eq. ( |T0D is equal to the third component of the vector density which is exactly zero for 
time reversal invariant wave functions. 

The above symmetry properties lead to extremely small overlaps between the four wave 
functions generated by the O^ symmetry operations: they are lower than 10 -5 for all the 
wave functions that we have used in our GCM calculation with a collective path defined 
along the z-axis. Our calculation has been performed without rotation. Two factors may 
increase these overlaps when the nucleus rotates. The first factor is the change in the spin 
content of the wave functions due to the Coriolis coupling. The second factor is the slight 
breaking of axial symmetry due to rotation. The importance of these effects remains to be 
studied quantitatively. 



IV. CONCLUSIONS 

In this work we have analyzed the presence of hexadecapole correlations in the mean 
fields created by Skyrme like effective interactions. Hexadecapole deformations have been 
introduced by the addition of hexadecapole constraints aligned either along the symmetry 
axis or the axis of rotation. 

The static effects of these correlations turn out to be very weak in SD 194 Hg and do 
not seem to be enhanced by rotation. We have also tried to put into evidence dynamical 
effects thanks the generator coordinate method. The GCM ground states obtained along 
the three collectives paths introduced do not show any fingerprint of polarisation due to 
the hexadecapole modes. The first GCM excited modes appear also rather high in energy, 
at 3-4 MeV. The lowest energy has been obtained in a GCM calculation along an x-path 
with relaxed quadrupole moment. We have interpreted this result as being due mainly 
to quadrupole correlations. However, it remains to be verified in a more complete study 
whether hexadecapole correlations are enhanced by a direct coupling between quadrupole 
and hexadecapole modes. 
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FIGURES 



FIG. 1. Potential energy curve as a function of the Q44 moment (top) and the Q44 moment 
(bottom). The solid line represents the result obtained from particle- number-projected wave func- 
tions without rotation and the dashed line represents the HFBLN results at 1=42. The dash-dotted 
line (bottom) represents the 1=0 results with a fixed quadrupole moment. 

FIG. 2. Mass quadrupole moment, (r 2 Y 2 o)> versus Q 44 in the nonrotating case (solid line) and 
at 1=42 (dashed line). 

FIG. 3. Squared moduli of the collective GCM wave functions plotted versus the Q| 4 moment 
(top) and the QJ 4 moment (bottom). Solid and dashed lines represent the results obtained with 
unconstrained and constrained quadrupole moments, respectively. 
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TABLES 



n 


Qaa 


Energy 


e* 


Qaa 


Energy 


e* 


Qaa 


Energy 


e* 


1 


3.5 


-1520.12 


0.00 


21.2 


-1520.15 


0.00 


20.4 


-1520.36 


0.00 


2 


-27.0 


-1516.10 


4.02 


20.1 


-1516.95 


3.21 


23.6 


-1516.17 


4.19 


3 


31.9 


-1514.47 


5.65 


15.5 


-1515.19 


4.97 


24.1 


-1515.14 


5.22 



TABLE I. The GCM results for the three collective paths specified by Q44, Q44 and Q44, where 
superscript c denotes the collective path with constrained quadrupole moment. Excitation energies 
are denoted by e*. Energies are given in MeV units and hexadecapole moments in fm 4 for the <5| 4 
path and in 1000 fm 4 in the two other cases. 



a 


X 


y 


z 


1 


+ 


+ 


Pk 


2 






Pk 


3 




+ 


-Pk 


4 


+ 




-Pk 



TABLE II. Parities of the components of a function of parity p^, with respect to the 
x = 0, y = and z = planes. 
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